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An  overview  of  GNSS  remote  sensing 
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Abstract 

The  Global  Navigation  Satellite  System  (GNSS)  signals  are  always  available,  globally,  and  the  signal  structures  are 
well  known,  except  for  those  dedicated  to  military  use.  They  also  have  some  distinctive  characteristics,  including 
the  use  of  L-band  frequencies,  which  are  particularly  suited  for  remote  sensing  purposes.  The  idea  of  using  GNSS 
signals  for  remote  sensing  -  the  atmosphere,  oceans  or  Earth  surface  -  was  first  proposed  more  than  two  decades 
ago.  Since  then,  GNSS  remote  sensing  has  been  intensively  investigated  in  terms  of  proof  of  concept  studies,  signal 
processing  methodologies,  theory  and  algorithm  development,  and  various  satellite-borne,  airborne  and  ground-based 
experiments.  It  has  been  demonstrated  that  GNSS  remote  sensing  can  be  used  as  an  alternative  passive  remote 
sensing  technology.  Space  agencies  such  as  NASA,  NOAA,  EUMETSAT  and  ESA  have  already  funded,  or  will  fund  in 
the  future,  a  number  of  projects/missions  which  focus  on  a  variety  of  GNSS  remote  sensing  applications.  It  is 
envisaged  that  GNSS  remote  sensing  can  be  either  exploited  to  perform  remote  sensing  tasks  on  an  independent 
basis  or  combined  with  other  techniques  to  address  more  complex  applications.  This  paper  provides  an  overview 
of  the  state  of  the  art  of  this  relatively  new  and,  in  some  respects,  underutilised  remote  sensing  technique.  Also 
addressed  are  relevant  challenging  issues  associated  with  GNSS  remote  sensing  services  and  the  performance 
enhancement  of  GNSS  remote  sensing  to  accurately  and  reliably  retrieve  a  range  of  geophysical  parameters. 


Review 

Recent  developments  in  both  ground-based  and  satellite- 
borne  geospatial  infrastructure  have  opened  exciting  new 
opportunities  for  geodesists  to  contribute  to  ‘big  issues' 
such  as  weather,  climate,  global  warming,  environment 
and  sustainability.  The  ability  to  remotely  sense  the  at¬ 
mosphere  using  geodetic  techniques  has  dramatically 
improved  over  the  past  decade,  primarily  as  a  result  of 
advances  in  satellite-borne  technologies,  large-scale  and 
dense  geodetic  Global  Navigation  Satellite  System  (GNSS) 
networks,  new  dedicated  space  missions  and  develop¬ 
ments  of  new  algorithms  and  innovative  methodologies. 
Atmospheric  sounding  is  a  new  area  of  GNSS  applications 
based  on  the  analysis  of  radio  signals  from  GNSS  satel¬ 
lites,  which  are  refracted  as  they  pass  through  the  atmos¬ 
phere  and  can  give  information  on  its  physical  properties. 
Depending  on  the  way  the  signals  (direct  or  reflected)  are 
used,  GNSS  remote  sensing  can  be  broadly  divided  into 
two  categories:  GNSS  reflectometry  (GNSS-R)  and  GNSS 
refractometry  (known  as  GNSS  radio  occultation  or 
GNSS-RO).  GNSS-R  involves  analysing  measurements  of 
GNSS  signals  reflected  from  the  Earths  surface,  while 
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GNSS-RO  utilises  measurements  of  GNSS  signals  refracted 
by  the  atmosphere  when  the  slanted  propagation  path  is 
close  to  the  Earths  limb.  In  the  following  sections  more 
details  about  the  theory  and  practice  of  GNSS  remote 
sensing  are  provided.  Interested  readers  are  also  referred 
to  [1,2]  for  more  background  information. 

Atmospheric  sensing  via  GNSS  meteorology 

The  use  of  GNSS  receivers  installed  on  low  Earth  orbit 
(LEO)  satellites  to  measure  the  refracted  radio  signals 
(GNSS-RO)  coupled  with  ground-based  atmospheric 
sounding  techniques  that  use  continuously  operating 
reference  station  (CORS)  networks  for  the  determin¬ 
ation  of  water  vapour  (WV)  variability  form  the  basis 
for  GNSS  meteorology.  The  dynamics  of  WV  have  a 
strong  influence  on  weather  and  climate  due  to  the  large 
energy  transfers  in  the  hydrological  processes.  This  is 
particularly  so  during  the  formation  and  life  cycle  of  se¬ 
vere  mesoscale  convective  storm  and  precipitation  sys¬ 
tems.  Contrary  to  its  importance,  WV  remains  poorly 
understood  and  inadequately  measured  both  spatially 
and  temporally,  especially  in  the  southern  hemisphere, 
where  meteorological  data  are  sparse. 
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GNSS  atmospheric  remote  sensing  using  LEO  satellites 
such  as  CHAllenging  Minisatellite  Payload  (CHAMP), 
Gravity  Recovery  and  Climate  Experiment  (GRACE)  and 
Constellation  Observing  System  for  Meteorology,  Iono¬ 
sphere,  and  Climate  (COSMIC)  is  a  robust  technique  for 
profiling  the  atmosphere  from  LEO  satellite  orbit  heights 
down  to  the  Earths  surface  [3-5].  Many  studies  have  dem¬ 
onstrated  that  the  GNSS-RO  technique  is  able  to  measure 
the  Earths  atmospheric  parameters  with  unprecedented 
high  accuracy,  high  resolution  and  global  coverage  [6-10]. 
RO  has  been  considered  as  a  new  observation  type  for 
meteorology  which  has  the  capability  to  provide  signifi¬ 
cant  information  on  the  thermodynamic  state  of  the  at¬ 
mosphere  and  to  allow  improvements  in  atmospheric 
analysis  and  prognosis.  GPS  RO  has  been  operationally 
used  in  numerical  weather  prediction  (NWP)  models  for 
weather  forecasting  by  major  weather  forecasting  centres 
such  as  those  in  the  USA,  Europe  and  Australia,  since 
2006.  The  beneficial  impact  of  such  data  assimilation  has 
been  quantified  in  several  studies  [11-14],  and  it  is  signifi¬ 
cant  particularly  in  data-void  regions  (e.g.  southern  hemi¬ 
sphere).  GNSS-RO  has  been  shown  to  significantly  reduce 
both  the  forecast  errors  in  the  ECMWF  system  and  NWP 
model  biases.  It  has  become  the  independent  sensor  tech¬ 
nique  with  the  fifth  highest  impact  in  operational  weather 
forecasting,  among  the  24  available  systems  [12-15].  Most 
recently,  Le  Marshall  et  al.  [11]  have  successfully  assimi¬ 
lated  MetOp,  GRACE  and  COSMIC  RO  data  into  the 
Australian  weather  forecasting  model  (ACCESS)  using 
4DVAR,  and  a  reliability  improvement  of  10  h  has  been 
achieved  in  a  1-  to  5 -day  period.  GPS  RO  has  been  offi¬ 
cially  used  in  operational  Australian  ACCESS  system  since 
late  2011. 

Building  on  the  success  of  the  COSMIC  RO  mission,  a 
follow-on  mission  (COSMIC-2)  was  recently  approved 
by  the  USA  and  Taiwan  that  will  launch  six  satellites 
into  low-inclination  orbits  in  late  2015  and  another  six 
satellites  into  high-inclination  orbits  in  early  2018.  Up  to 
12,000  high-quality  profiles  per  day  will  be  tracked  when 
it  is  fully  deployed.  This  revolutionary  increase  in  the 
number  of  atmospheric  and  ionospheric  observations 
will  greatly  benefit  the  research  and  operational  communi¬ 
ties  and  will  open  up  unprecedented  research  opportun¬ 
ities  for  hurricane  analysis  and  prediction,  and  studies  of 
Madden- Julian  Oscillation  and  climate  processes,  to  name 
just  a  few.  In  addition,  algorithm  development,  optimisa¬ 
tion  and  refinement,  error  reduction,  climate  analysis  and 
benchmarking,  and  polar  applications  are  all  frontier  areas 
of  RO  research  [15-21]. 

Vertical  gradients  of  the  atmospheric  properties  are  re¬ 
trieved  from  precise  measurements  of  the  radio  propaga¬ 
tion  path  between  a  LEO-based  GPS  receiver  and  an 
occulting  GPS  satellite  transmitter  pair.  A  radio  occult- 
ation  event  (ROE)  occurs  when  a  GPS  satellite  is  setting 


or  rising  behind  the  Earths  limb  as  viewed  by  a  receiver 
onboard  a  LEO  satellite  as  shown  in  Figure  1.  During  a 
ROE,  the  transmitted  signals  from  the  GPS  to  the  LEO 
satellite  are  delayed.  This  delay  can  be  converted  into  a 
profile  of  bending  angles  and  consequently  processed  to 
produce  a  refractivity  index  profile.  The  refractivity  is  a 
function  of  the  electron  density  in  the  ionosphere,  and 
temperature,  pressure  and  WV  in  the  atmosphere.  The 
fundamental  GPS  observation  equation  based  upon  which 
atmospheric  parameters  are  retrieved  can  be  expressed  as 
follows: 

L  =  p  +  c(Stc\^-StTe\)—c{S  Tc\fc-5Tre\)  +  ALA 

+  excess  +  AT?  +  £  (l) 

where  L  is  the  raw  GPS  carrier  phase  measurement,  p 
is  the  geometric  distance  between  the  LEO  and  GPS 
satellites,  c  is  the  speed  of  light,  (8tc n<,  8tre\)  are  the  er¬ 
rors  related  to  the  clock  and  relativistic  effects  of  LEO, 
(8Tclk,  8Trel)  are  the  errors  related  to  GPS  satellite 
clock  and  relativistic  effects,  N  is  the  integer  ambigu¬ 
ity,  A  is  the  wavelength  of  the  GPS  L-band  signal, 
A^excess  is  the  excess  phase  along  the  radio  signal 
propagation  path  which  is  the  phase  increment  due  to 
the  signal  propagation  delay,  A R  is  the  general  relativistic 
effect,  and  e  is  the  measurement  noise.  Figure  2  shows  the 
flowchart  of  a  typical  GNSS-RO  atmospheric  retrieval 
process,  and  Figure  3  shows  a  typical  1-day  RO  event  distri¬ 
bution  across  the  globe  under  a  full  COSMIC  constellation. 

The  L-band  radio  signals  broadcast  by  the  GNSS  satel¬ 
lites  are  affected  by  both  the  ionospheric  and  tropospheric 
effects  as  they  travel  through  the  atmosphere  to  the  per¬ 
manent  ground-based  GNSS  receivers.  Using  a  dense 
CORS  network,  the  impact  of  the  neutral  atmosphere,  i.e. 
tropospheric  delays,  can  be  estimated  as  a  by-product  of 
the  geodetic  data  processing.  This  ground-based  GNSS 
meteorology  is  currently  used  to  produce  two  different 
data  products.  Firstly,  the  integrated  precipitable  WV  in 
the  zenith  direction  above  each  GPS  station  can  be  accur¬ 
ately  estimated  using  path  delay  measurements  between  a 
receiver  and  GPS  satellites.  Secondly,  the  spatial  and  tem¬ 
poral  distribution  of  wet  refractivity  can  be  reconstructed 
using  accurate  measurements  of  path  delays  between  each 
satellite-and-receiver  link  over  the  GNSS  network.  This 
method  is  called  GNSS  tomography  which  is  an  emerging 
method  for  4D  tropospheric  reconstruction  using  a  local 
GNSS  CORS  network.  Integrated  measurements  of  wet 
delays  are  extracted  from  each  satellite- to- station  link, 
which  are  used  as  the  primary  observable  through  a  finite 
3D  grid  field  to  estimate  the  time-dependent  3D  structure 
of  WV.  An  integrated  approach  is  usually  used  in  the 
tomographic  reconstruction  process  by  incorporating 
multiple  sources  of  measurements  from  GNSS  CORS, 
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Figure  1  Illustration  of  the  geometry  of  the  direct,  reflected  and  refracted  GNSS  signals  and  their  use  for  different  applications. 


automated  weather  stations,  radiosonde,  RO  and  NWP 
models. 

A  variety  of  research  has  been  conducted  based  on  both 
real  and  simulated  data.  For  example,  [22]  presented  a 
simulation  and  real  data  investigation  into  4D  tropo¬ 
spheric  tomography  using  GPS  observations  and  a  sto¬ 
chastic  approach.  The  real  data  analysis  used  the  Kilauea 
network  in  Hawaii  with  the  voxel  grid  discretisation  cover¬ 
ing  an  area  of  400  km2  and  15  km  in  height.  A  Kalman  fil¬ 
tering  technique  was  used  for  estimating  the  evolution  of 
dynamically  changing  parameters.  An  alternate  method 
presented  in  [23]  used  raw  GPS  phase  measurements 


Figure  2  Typical  data  processing  steps  involved  in  GNSS-RO 
atmospheric  retrievals. 


directly  for  the  reconstruction  of  WV  distribution.  This 
method  was  tested  using  simulations  where  the  horizontal 
and  vertical  voxel  resolution,  observation  noise  and  add¬ 
itional  GNSS  constellations  were  varied  and  optimised. 
Results  indicated  that  high  accuracy  can  be  achieved  in 
the  modelling  of  wet  refractivity  in  lower  layers.  However, 
its  limitation  was  in  resolving  layers  above  4,000  m 
altitude.  Nilsson  and  Gradinarsky  [24]  investigated  the 
benefits  of  GPS  tomography  for  modelling  mesoscale 
precipitation  systems.  The  model  was  designed  on  the 
premise  that  the  grid  system  would  move  in  relation  to 
the  vector  of  horizontal  wind  velocity.  It  was  concluded 
that  GPS  tomographic  solutions  can  improve  the  initial 
WV  field  for  short-term  NWP  predictions  and  forecasting. 
Bender  and  Raabe  [25]  extensively  analysed  the  precondi¬ 
tions  and  constraints  of  ground-based  GPS  WV  tomog¬ 
raphy  and  suggested  that  the  spatial  resolution  is  limited 
by  the  density  of  the  ground-based  GPS  networks  and  the 
number  of  available  GNSS  satellites.  To  optimise  horizon¬ 
tal  resolution,  the  horizontal  size  of  the  humidity  field 
voxels  would  be  restricted  to  the  mean  interstation  dis¬ 
tance  of  the  GPS  network.  Another  important  parameter 
presented  is  the  elevation  cutoff  angle;  the  number  of 
available  rays  increases  heavily  with  a  decreasing  elevation 
cutoff  angle  and  in  return  provides  much  more  informa¬ 
tion  on  the  lower  troposphere.  Rohm  and  colleagues 
[26,27]  presented  investigations  into  the  concept  of  near- 
real-time  atmospheric  model  construction  based  on  the 
GNSS  and  meteorological  observations.  This  nowcasting 
approach  provides  accurate  comparisons  with  NWP 
analysis  models.  Yuan  et  al.  [28]  recently  investigated 
real-time  retrieval  using  precise  point  positioning  tech¬ 
nique  and  good  results  were  obtained.  Choy  et  al.  [29] 
investigated  the  propagation  and  distribution  of  WV 
during  the  Melbourne  severe  weather  event  in  2010, 
while  Manning  and  colleagues  [30,31]  investigated  4D 
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Figure  3  One-day  RO  event  distribution  across  the  globe  under  a  full  COSMIC  constellation.  GPS  RO  sounding  locations  for  COSMIC 
(in  green)  and  radiosonde  sites  (in  red)  in  a  period  of  24  h  (after  UCAR). 


Occultation  Locations  for  COSMIC,  6  S/C,  6  Planes,  24  Hrs 


methods  of  reconstructing  wet  refractivity  fields  using 
GPS  tomography  for  the  Australian  region.  A  comprehen¬ 
sive  investigation  was  carried  out  in  terms  of  sensitivity 
analysis  of  multiple  measurements,  voxel  optimisation, 
parameterisation  method  and  Kalman  filtering  using  two 
severe  storm  events  in  Melbourne  as  case  studies.  Bennitt 
and  Jupp  [32]  studied  the  data  assimilation  of  GPS  zenith 
total  delays  into  the  Met  Office  NWP  models,  and  a  slight 
positive  overall  impact  and  significant  improvements  to 
short-range  forecasts  in  some  cases  were  demonstrated. 

To  summarise,  GNSS  is  now  an  established  atmos¬ 
pheric  observing  system  which  can  accurately  sense  WV, 
the  most  abundant  greenhouse  gas,  accounting  for  60  to 
70%  of  atmospheric  warming.  Severe  weather  forecasting 
is  challenging,  in  part  due  to  the  high  temporal  and 
spatial  variation  of  atmospheric  WV.  WV  content  is 
undersampled  in  the  current  meteorological  and  climate 
observing  systems;  obtaining  and  exploiting  more  high- 
quality  humidity  observations  is  essential  to  NWP  data 
assimilation,  nowcasting  and  climate  monitoring.  It  is 
expected  that  the  GNSS  meteorology  technology  will 
significantly  advance  our  knowledge  of  Earths  atmospheric 
structure  and  processes,  including  severe  weather  phe¬ 
nomena  and  climate  change,  in  particular  in  data-sparse 
areas  such  as  the  southern  hemisphere,  and  oceanic  and 
polar  regions. 

Ocean  observation 

The  complex  temporal  and  spatial  variations  which  shape 
the  ocean  surface  are  amenable  to  study  using  remote 
sensing  methods  because  the  sea  surface  both  emits  and 
strongly  reflects  electromagnetic  radiation.  Static  ocean 
surface  undulations  are  dominated  by  Earths  geoid  and  by 
its  topographic  response  to  the  mean  ocean  circulation, 
while  temporal  departures  from  this  surface  are  forced  by 
weather  and  climate-related  changes  in  ocean  temperature 


and  salinity,  and  surface  wind  stress.  Lunar  and  solar 
gravitational  tides,  and  their  harmonics,  contribute 
additional  variability,  on  daily  and  longer  timescales. 
These  sea  level  changes,  which  also  drive  ocean  cur¬ 
rents,  have  been  sensed  remotely  using  nadir-viewing 
satellite-borne  radar  altimeters,  almost  continuously 
over  the  last  four  decades.  Observations  from  these  altim¬ 
eters  provide  vital  information  for  studies  of  deep-ocean 
circulation  and  boundary  currents,  the  mid-ocean  gyres, 
tsunamis  and  ocean  currents  on  synoptic  to  global  scales, 
while  mesoscale  (approximately  100  km)  ocean  eddies  are 
only  marginally  resolved. 

Finer-scale  topographic  changes  associated  with  ocean 
surface  roughness,  or  sea  state’,  are  dominated  by  locally 
forced  wind-waves  and  swell  from  distant  storms  occur¬ 
ring  on  time  and  space  scales  exceeding  a  few  hours  and 
kilometres.  In  addition  to  their  effects  on  shipping  and 
other  human  activities,  these  contribute  to  ocean  mix¬ 
ing,  which  crucially  affects  the  exchange  of  momentum, 
heat,  salt  and  various  gases,  such  as  C02,  across  the  air- 
sea  interface,  with  important  implications  for  weather 
and  climate  studies.  Sea  state  also  affects  the  interaction 
of  electromagnetic  radiation  with  the  sea  surface  through 
its  impact  on  the  effective  emissivity  and  reflectivity  of  the 
surface.  By  using  suitable  remote  sensing  methods,  its  ef¬ 
fects  can  be  used  to  deduce  sea  state  or  mean  square 
slope.  Coupled  with  a  wind-wave  spectral  model,  the  sur¬ 
face  wind  field  can  also  be  inferred.  Alternatively,  these  ef¬ 
fects  can  be  used  to  derive  corrections  for  roughness 
influence  on  retrievals  of  other  geophysical  parameters 
from  microwave  radiometers,  such  as  sea  surface  salinity. 
Since  sea  state  modifies  the  waveform  of  the  radar  reflec¬ 
tion,  nadir-viewing  (monostatic)  radar  altimeters  are  also 
used  to  estimate  wave  height  and  wind  speed,  with  high 
spatial  resolution  (approximately  7  km)  along-track  but 
poorer  resolution  (approximately  300  km)  across-track 
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and  with  an  approximately  10-day  repeat  cycle  [33].  How¬ 
ever,  over  the  last  three  decades,  a  new  approach  to  meas¬ 
uring  sea  surface  height,  wind  speed  and  sea  state  bistatic 
radar,  offering  higher  spatial  and  temporal  resolution,  and 
utilising  GNSS  satellites  as  transmitters  has  been  proposed 
and  developed. 

Ocean  observation  using  GNSS  signals  was  first  pro¬ 
posed  as  the  Passive  Reflectometry  and  Interferometry 
System  (PARIS)  concept  to  perform  mesoscale  ocean 
altimetry  [34],  and  an  in-orbit  demonstrator  has  been 
designed  to  implement  this  concept  [35].  Recent  develop¬ 
ments  in  GNSS-R  altimetry  include  experimental  ground- 
based  and  airborne  campaigns,  performance  studies, 
refinements  in  error  budget  estimates  and  alternative 
signal  processing  methods  for  retrieving  sea  surface 
height  estimates  from  GNSS-R  instrumentation.  The  new 
methods  exploit  a  variety  of  techniques  for  enhancing  per¬ 
formance  characterised  by  improvements  to  spatial  and 
temporal  resolutions,  signal  to  noise  ratio  and/or  accuracy 
and  precision.  With  respect  to  error  characterisation, 
Martin-Neira  et  al.  [36]  studied  the  effects  of  errors  in 
delay  drift  compensation  on  the  waveform  and  its  impact 
on  the  altimetric  accuracy  and  precision,  while  Pascual 
et  al.  [37]  used  simulated  GPS  and  Galileo  signals  to  inves¬ 
tigate  the  impact  of  bandwidth  and  low-power  transmit¬ 
ters  on  height  accuracy  and  precision.  Rius  et  al.  [38] 
investigated  the  scatterometric  and  specular  delay  observ¬ 
ables  and  considered  various  sources  of  uncertainty  in¬ 
cluding  observation  statistics  and  model  parameterisation 
or  calibration  error.  Methodology  studies  include  an  in¬ 
vestigation  of  model-free  retrieval  of  sea  surface  height 
using  a  power  ratio  method  combined  with  cost  function 
optimisation  when  the  sea  surface  is  rather  rough  [39].  A 
semi-codeless  processing  approach  [40]  showed  how  to 
determine  the  phase  sign  variation  of  encrypted  P-code 
GNSS  signals  and  combine  it  with  the  CA  signal  stream 
to  enhance  signal  to  noise  ratio  in  comparison  with  that 
achieved  by  the  interferometric  method.  An  investigation 
of  Doppler  multilook  processing  [41]  realised  enhanced 
spatial  and  temporal  resolution  by  intra-  and  inter- 
Doppler  tracking  and  averaging  of  delay  waveforms.  The 
method  demonstrated  25  to  30%  improvement  in  alti¬ 
metric  precision  over  the  interferometric  method  and 
increased  flexibility  in  trading  off  spatial  resolution  and 
height  estimation  accuracy.  In  a  ground-based  study  re¬ 
ported  in  2001  [33],  retrieval  of  Crater  Lake  surface 
height  to  2-cm  precision  was  achieved  using  GPS  altim¬ 
etry  from  relatively  low  satellite  elevations.  A  year  later, 
5-cm  precision  GPS  ocean  altimetry  was  achieved  from 
an  aircraft  [42].  A  ground-based  experiment  implement¬ 
ing  the  PARIS  concept  in  the  form  of  an  X-band  instru¬ 
ment  [43]  takes  advantage  of  high  signal  to  noise  ratio 
(SNR)  associated  with  X-band  signals  of  opportunity 
from  digital  satellite  TV  broadcasts.  Useful  overviews 


including  summaries  of  these  and/or  other  GNSS  altim¬ 
etry  studies  appear  in  [2,44,45],  along  with  an  assess¬ 
ment  of  the  overall  performance  of  GNSS-R  altimetry 
techniques. 

Using  airborne  data  collected  in  1997,  Garrison  et  al. 
[46]  first  demonstrated  the  potential  for  retrieving  sea 
state  from  the  dependence  of  the  GNSS  sea  surface  reflec¬ 
tion  waveform  width  on  wind  speed.  This  waveform  rep¬ 
resented  the  cross-correlation  between  a  replica  of  the 
direct  signal  and  the  signal  reflected  off  the  sea  surface,  as 
a  function  of  delay.  Coupled  with  knowledge  of  the  GPS 
navigation  code  autocorrelation  function  (which  is  usually 
triangular  in  shape),  this  gives  an  estimate  of  the  radar 
cross  section  that  is  closely  related  to  the  ocean  surface 
mean  square  slope  and  hence  surface  roughness.  During 
the  period  2003  to  2011,  the  UK-Disaster  Monitoring 
Constellation  (UK-DMC)  satellite  carried  an  experimental 
GNSS  reflectometer,  which  yielded  sea  state  and  wind 
speed  estimates  for  representative  cases  of  low  to  moder¬ 
ate  wind  strength  [47].  While  suitable  for  scatterometric 
measurements,  this  instrument  was  not  optimised  for  alti¬ 
metric  measurements. 

Since  Garrison  s  first  experiment,  a  variety  of  techniques 
have  been  developed  for  retrieving  surface  mean  square 
slope  and  wind  speed,  and  many  have  been  tested  using 
ground  and  airborne  field  campaigns.  The  following  dis¬ 
cussion  of  some  key  studies  exemplifies  the  progress  made 
but  is  not  meant  to  be  comprehensive:  Retrieval  of  ocean 
wind  speeds  from  a  balloon  in  the  stratosphere  [48]  pro¬ 
duced  mean  square  slope  and  corresponding  wind  speed 
estimates  departing  from  radiosondes  and  satellite  altime¬ 
ters  and  a  scatterometer  by  at  most  1  m/s  in  relatively 
calm  winds  (less  than  about  9  m/s).  Retrieval  of  hurricane 
winds  using  GPS  reflectometry  was  first  demonstrated  by 
Katzberg  et  al.  [49]  who  later  calibrated  retrieved  hurri¬ 
cane  winds  using  dropsondes  [50,51].  The  GNSS-R  re¬ 
flectometer  was  found  to  underestimate  observed  winds, 
but  a  standard  deviation  around  a  linear  regression  fit  to 
observations  of  about  5  m/s  was  found,  for  winds  ranging 
from  10  to  40  m/s.  Recently,  wind  speeds  have  been  re¬ 
trieved  using  GNSS-R  instrumentation  from  an  aircraft 
using  scattering  model-free  linear  regression  of  several 
Delay  Doppler  Map  (DDM)  observables  [52].  The  results 
obtained  using  satellite  triplets,  and  validated  using  coinci¬ 
dent  dropsonde  measurements,  produced  errors  of  up  to 
1  to  4  m/s  for  winds  of  5  to  25  m/s,  with  overall  accuracy 
and  precision  of  the  order  of  1  to  2  m/s. 

Methods  to  extract  directional  roughness  characteristics 
have  been  developed  and  tested  using  airborne  campaign 
data  [53-55].  Wind  direction  retrieval  is  performed  by 
accounting  for  differences  in  upwind/crosswind  surface 
slope  probability  distribution,  which  is  represented  by  an 
observable  asymmetry  in  the  glistening  zone.  The  particu¬ 
lar  technique  devised  by  Cardellach  and  Rius  [55]  which 
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does  not  require  an  assigned  probability  distribution 
function  can  resolve  the  180°  directional  ambiguity  in¬ 
herent  in  other  GNSS-R  and  radar  scatterometry  pro¬ 
cessing  methods.  Retrieval  of  roughness  characteristics 
using  relatively  powerful  signals  of  opportunity  from 
S-band  satellite  radio  broadcasts  produced  SNR  values 
greatly  exceeding  those  of  GPS  transmissions  [56]. 
These  broadcasts  can  provide  high-quality  GNSS-R  + 
data  along  the  coastal  margins  but  are  not  available 
over  the  deep  ocean. 

Useful  reviews  of  ground-based  and  airborne  GNSS-R 
campaigns,  geophysical  models  and  signal  processing 
techniques  for  a  variety  of  oceanographic  applications  are 
available  [2,44].  Campaigns  performed  over  ocean  surfaces 
have  spanned  locations  in  the  North  Sea,  the  Baltic  Sea, 
the  Mediterranean  Sea  and  the  Atlantic  Ocean  [44],  and 
the  Tasman  Sea  in  the  western  South  Pacific  [39].  Newly 
developed  GNSS-R  signal  processing  techniques,  exploiting 
favourable  sampling  geometries  or  instrument  refinements, 
can  now  be  used  to  extract  such  factors  as  significant  wave 
height  and  roughness  field  directionality  [38,54]  or  to  re¬ 
late  ocean  roughness  to  surface  brightness  temperature 
measured  at  microwave  frequencies  [57,58].  The  possibil¬ 
ity  of  retrieving  ocean  surface  temperature  and  salinity 
from  permittivity  estimates  derived  from  polarimetric 
GNSS-R  measurements  has  also  been  investigated  [59]. 
Such  a  development  could  complement  the  retrieval  of 
sea  surface  salinity  from  current  L-band  radiometry  mis¬ 
sions  (ESAs  SMOS  and  NASA's  Aquarius),  which  is 
strongly  affected  by  the  influence  of  roughness  effects  on 
emissivity. 

A  well-established  approach  to  estimate  sea  state  and 
wind  speed  using  GNSS-R  measurements  is  based  on 
the  theoretical  model  of  [60],  which  solves  the  bistatic 
radar  equation  of  the  time-delayed  scattered  power  for 
the  radar  cross  section  a0: 


(\nr,fc)\2) 


G2{p)\S(fc-fD(fi))\2A2 

T2cr  (t-{Rq(?)-R(P))/c) 
a  4  nR2(p)R2(p) 


ff0(p)  d2p 


(2) 


Here,  t  is  the  time  delay  relative  to  the  specular  reflec¬ 
tion;  Tt  is  the  averaging  time;  A  is  the  reflecting  surface 
area;  p  is  the  reflecting  surface  position  vector;  and  R0 
and  R  are  the  distances  from  a  surface  point  to  the 
transmitter  and  receiver,  respectively;  and  fc  and  fD  are 
the  frequencies  of  the  replica  signal  and  the  reflected 
Doppler  signal,  respectively.  The  terms  G,  S  and  A  rep¬ 
resent  the  antenna  pattern,  the  sine  function  for  the 
Doppler  mismatch  and  the  navigation  signal  correlation 
function,  respectively.  In  general,  this  equation  is  solved 
in  2D  delay  Doppler  space  using  various  geometric  de¬ 
scriptions  of  the  corresponding  DDM  [61].  However,  it 
is  often  simpler  and  faster  to  perform  retrievals  in  ID 


delay  space  by  deriving  the  corresponding  delay  wave¬ 
forms.  In  the  limit  set  by  the  geometric  optics  theory,  o0 
is  proportional  to  the  probability  density  function  of 
surface  gravity  wave  slopes  from  which  the  mean  square 
slope  (MSS)  may  be  derived.  A  semi-empirical  wind-wave 
spectrum  [62]  is  typically  used  to  infer  wind  speed  at 
standard  meteorological  height  from  the  derived  MSS. 

An  example  of  wind  speed  estimation,  based  on  the 
model  fitting  method,  appears  in  Figure  4  [63].  Five  the¬ 
oretical  ID  delay  waveforms  are  shown  corresponding  to 
five  different  wind  speeds  (3,  4,  5,  6  and  7  m/s),  with  a 
measured  delay  waveform  from  GPS  satellite  PRN#13 
superimposed.  The  data  were  collected  by  an  airborne 
receiver  at  an  altitude  of  about  3  km.  The  measured 
waveform  is  the  result  of  coherent  integration  of  1-ms 
intermediate  frequency  (IF)  signals  followed  by  inco¬ 
herent  integration  of  1,000  of  these  1-ms  integrals.  The 
measured  waveform  matches  well  with  the  theoretical 
waveform  of  wind  speed  4  m/s,  which  closely  matches 
the  observed  wind  speed  of  4.4  m/s. 

Recent  successes  in  global  mapping  of  sea  surface  sal¬ 
inity  using  L-band  radiometers  aboard  ESAs  SMOS  and 
NASA's  Aquarius  satellites  provide  additional  motivation 
for  developing  techniques  to  deduce  corrections  for  rough¬ 
ness  influence  on  L-band  emissivity.  Several  airborne  mis¬ 
sions  have  demonstrated  this  potential,  and  the  planned 
launch  of  NASA's  Cyclone  Global  Navigation  Satellite 
System  (CYGNSS)  L-band  reflectometer  constellation  in 
2016,  with  the  goal  of  mapping  surface  winds  at  a  high 
temporal  resolution  inside  tropical  cyclones,  will  allow  this 
concept  to  be  tested  in  space.  Recent  extensions  of  GNSS- 
R  to  communications  satellites  operating  on  higher  fre¬ 
quency  bands,  with  improved  signal  strength,  demonstrate 
enhanced  precision  of  ocean  reflectometry  measurements 
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with  the  possibility  of  expanded  coverage  in  coastal  re¬ 
gions  [43,56,64]. 

Land  applications 

Knowledge  of  soil  moisture  content  is  critical  for  drought 
and  irrigation  management,  so  as  to  increase  crop  yields 
and  to  gain  a  better  understanding  of  natural  processes 
linked  to  the  water,  energy  and  carbon  cycles.  Soil  mois¬ 
ture  can  be  estimated  by  analysing  surface-scattered  sig¬ 
nals  transmitted  and  received  by  active  radar  sensors  or 
natural  surface  emission  detected  by  microwave  radiome¬ 
ters.  Alternatively,  GNSS-R  signals  can  be  employed  to  es¬ 
timate  soil  moisture.  Since  the  soil  dielectric  constant  and 
reflectivity  depend  on  soil  moisture,  the  variation  in  the 
reflected  GNSS  signal  power,  or  SNR,  indicates  changes  in 
soil  moisture  [65].  If  feasible,  the  empirical  dielectric  model 
based  method  in  [66]  may  be  applied  to  GNSS -based  soil 
moisture  retrieval.  In  this  method,  the  transition  moisture 
parameter  ( 6t )  is  defined  as 

6>t  =  0.165  +  0.49  flwt  (3) 

where  O^t  is  the  wilting  point  moisture  that  is  a  function 
of  the  percentage  of  the  clay  and  sand  contents.  Depend¬ 
ing  on  whether  the  soil  moisture  is  greater  or  less  than 
the  transition  moisture,  a  specific  formula  can  be  used 
to  calculate  the  soil  moisture.  This  model  employs  the 
mixing  of  either  the  dielectric  constants  or  the  refraction 
indices  of  ice,  water,  rock  and  air  and  treats  the  transi¬ 
tion  moisture  value  as  an  adjustable  parameter.  The  di¬ 
electric  constant  of  the  soil  can  be  estimated  by 
measuring  the  surface  reflectivity.  The  model  is  quite 
general  since  it  was  developed  by  considering  many 
different  types  of  soils. 

An  alternative  GNSS-R  technique  for  soil  moisture  esti¬ 
mation  is  the  interferometric  method  in  which  the  power 
associated  with  the  coherent  sum  of  the  direct  and 
reflected  GNSS  signals  is  exploited.  After  continuously 
measuring  the  total  power  for  a  few  hours,  the  location  of 
the  power  notch  over  satellite  elevation  can  be  observed, 
and  the  soil  moisture  may  be  estimated  [67].  This  notch- 
based  approach  is  limited  to  the  case  where  a  vertically 
polarised  antenna  is  used  and  fixed  as  horizon-looking.  It 
would  not  be  applicable  to  standard  geodetic  antennas.  It 
is  an  advantage  to  use  existing  GPS  receivers,  installed  pri¬ 
marily  for  geophysical  and  geodetic  applications,  to  esti¬ 
mate  soil  moisture  [66].  The  zenith-looking  antenna  at  a 
geodetic  reference  receiver  station  captures  both  direct 
and  reflected  GPS  signals,  although  the  reflected  signals 
are  received  at  a  negative  elevation  angle  where  the  an¬ 
tenna  radiation  gain  is  non-zero.  The  amplitude  (or  SNR) 
and  phase  (or  frequency)  of  the  multipath  signal  can  be 
used  to  infer  soil  moisture  information  [68].  These  re¬ 
ceivers,  if  exploited  effectively,  could  provide  a  global 


network  for  services  on  soil  moisture,  vegetation  and  snow 
cover  [69,70].  GNSS-based  soil  moisture  estimation  is  com¬ 
plicated  by  a  number  of  issues  including  surface  roughness, 
vegetation  canopy  and  variation  in  percentage  of  individual 
soil  components.  To  achieve  reliable  soil  moisture  estima¬ 
tion,  these  issues  must  be  taken  into  account  through  pro¬ 
cesses  such  as  modelling  and  compensation. 

Forest  change  detection  is  another  possible  application 
of  GNSS  remote  sensing.  Forest  change  can  provide  use¬ 
ful  information  about  global  climate  change  impacts  and 
carbon  storage,  and  knowledge  of  forest  change  is  vital 
for  effective  forest  management.  Received  signal  strength 
of  the  reflected  GNSS  signals  can  be  employed  to  distin¬ 
guish  forest  conditions  from  each  other,  since  different 
surface  covers  have  different  reflectivities.  A  number  of 
signal  strength  ranges  may  be  defined  to  be  associated 
with  a  group  of  surfaces  such  as  lake/river  water,  typical 
dense  forest  and  cleared  area  due  to  logging  [71].  Figure  5 
shows  four  ground  specular  reflection  tracks  of  the 
reflected  signals  from  four  GNSS  satellites  received  by  a 
receiver  on  an  aircraft.  The  tracks  are  colourised  by  the 
reflected  signal  power.  The  low  signal  power  corresponds 
to  the  dense  forest  areas,  while  the  high  signal  power  oc¬ 
curs  over  three  areas  marked  by  Al  (cleared  area),  A2 
(partially  cleared  area)  and  A3  (lake  water).  Accordingly,  a 
surface  can  be  classified  by  determining  within  which 
range  its  measured  strength  falls.  In  addition  to  signal 
strength,  other  signal  characteristics  such  as  those  related 
to  the  observed  correlation  waveform  may  be  used  to  en¬ 
hance  surface  change  detection  or  surface  classification. 
The  major  issue  related  to  surface  change  detection  is  the 
dependence  of  reflectivity  on  several  factors  such  as  soil 
moisture  and  surface  roughness.  Forest  change  may  be 
quantified  by  evaluating  forest  biomass,  which  depends  on 
the  volume  of  both  living  and  dead  trees  (leaves,  branches 
and  trunks).  The  received  signal  power  can  be  used  to  de¬ 
termine  the  scattering  coefficient  which,  based  on  simula¬ 
tion  [72],  should  be  a  function  of  the  biomass.  However,  it 
is  a  challenging  problem  to  derive  a  formula  that  accur¬ 
ately  describes  the  relationship  between  the  scattering  co¬ 
efficient  and  the  biomass. 

Cryosphere  applications 

GNSS-R  has  been  used  to  sense  sea  ice  and  snow,  mid- 
to  high-latitude  continental  snow  and  deep  sub-surface 
layers  of  Antarctica's  dry  snow.  Three  main  parameters 
of  sea  ice  (i.e.  thickness,  surface  roughness  and  ice  per¬ 
mittivity)  can  be  retrieved  with  GNSS-R.  These  parame¬ 
ters  can  be  combined  to  help  characterise  different  ice 
types  including  new  ice,  young  ice,  thin  first-year  ice, 
first-year  ice,  and  multiyear  ice.  Sea  ice  thickness  is  a  key 
parameter  for  classification  and  characterisation  of  sea 
ice  masses,  which  influence  the  temperature  and  circula¬ 
tion  pattern  of  both  the  ocean  and  atmosphere  and  thus 


Yu  et  al.  EURASIP  Journal  on  Advances  in  Signal  Processing  2014,  2014:134 
http://asp.eurasipjournals.com/content/201 4/1  /1 34 


Page  8  of  14 


can  be  used  for  analyses  of  the  Earths  climate.  In  the 
case  of  reflections  over  a  sea  ice  surface,  it  is  possible  to 
track  phase  delay  precisely  so  as  to  retrieve  the  ice  thick¬ 
ness  through  phase  altimetry.  ESA  initiated  a  GPS  sea-ice 
campaign,  and  a  ground-based  experiment  was  conducted 
in  Greenland  over  a  7-month  period  in  order  that  the 
complete  process  of  sea  ice  formation  and  melting  was 
monitored  [73,74].  The  time  series  of  the  retrieved  heights 
is  in  agreement  with  an  Arctic  tide  model,  and  there  is  a 
consistency  between  results  obtained  using  different  an¬ 
tenna  polarisations. 

Similar  to  sea  surface  wind  speed  retrieval,  model  fit¬ 
ting  between  the  theoretical  and  measured  delay  wave¬ 
forms  can  be  used  to  determine  the  roughness  of  an  ice 
surface  [75].  In  the  transition  from  newly  formed  sea  ice 
to  older  and  thicker  ice,  a  gradual  decrease  was  observed 
in  bulk  salinity  and  an  increase  in  surface  roughness. 
Roughness  may  also  be  inferred  by  investigating  the 
characteristics  of  the  phase  or  coherence  of  the  inter¬ 
ferometric  signal  [73,76].  Permittivity  of  sea  ice  is  an¬ 
other  important  parameter  which  can  be  estimated  by 
determining  the  ratio  between  the  co-polar  and  cross- 
polar  components  of  the  reflected  signal  [44,76]  or  by 
model  fitting  of  delay  waveform  [77]. 

Snow  is  a  critical  storage  component  in  the  hydrologic 
cycle  and  an  important  component  of  the  climate  sys¬ 
tem;  however,  current  in  situ  observations  of  snow  dis¬ 
tribution  are  sparse.  By  making  use  of  the  networks  of 
GPS  stations  established  for  geodetic  applications,  it  is 
possible  to  monitor  snow  distribution  on  a  global  scale. 
The  SNR  and  phase  of  the  multipath  signal  are  both  a 
function  of  the  antenna  height  above  the  ground.  The 
presence  of  a  snow  pack  around  the  antenna  reduces  the 


antenna  height  and  hence  changes  the  multipath  signal 
SNR  [78,79].  The  multipath  signal  SNR  is  a  periodic 
function  with  a  carrier  phase  given  by 

4tz7z  . 

<Pm=  ~Y~  Sm0  (4) 

where  h  is  the  antenna  height  relative  to  the  ground  sur¬ 
face,  X  is  the  wavelength,  and  6  is  the  satellite  elevation 
angle.  Differentiating  the  phase  with  respect  to  sin  6  will 
produce  the  signal  frequency 


As  a  consequence,  the  antenna  height,  and  hence  the 
snow  depth,  can  be  estimated  by  measuring  the  fre¬ 
quency  of  the  multipath  signal  SNR.  It  is  worth  men¬ 
tioning  that  GNSS  signal  reflection  off  Antarctic  dry 
snow  can  be  significantly  different  from  that  of  typical 
wet  snow  surfaces.  The  L-band  GNSS  signals  can  pene¬ 
trate  into  the  dry  snow  down  to  200  to  300  m,  produ¬ 
cing  a  complex  interference  pattern  induced  by  the 
multiple  reflections  occurring  at  different  layer  interfaces 
of  the  sheet  [80].  Radio-holographic  techniques  are  used 
on  each  lag  of  the  delay  waveform  to  identify  the  spec¬ 
tral  content  of  the  signal  and  to  identify  each  frequency 
component  to  different  snow  depths. 

Missions 

Useful  reviews  of  radio  occultation  missions  can  be  found 
in  [81]  and  [82].  The  concept  was  first  proved  by  the 
GPS/MET  (GPS/Meteorology)  satellite  mission  in  1995  to 
1997.  That  satellite  took  few  measurements  but  was  followed 
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by  the  more  productive  CHAllenging  Minisatellite  Payload 
(CHAMP)  and  Satellite  de  Aplicaciones  Cientificas-C 
(SAC-C)  satellites.  CHAMP  provided  8  years  of  radio  oc- 
cultation  data  consisting  of  around  440,000  measurements 
from  February  2001  to  October  2008  [83].  Those  missions 
led  to  the  launch  in  2006  of  a  six-satellite  constellation 
FORMOSAT-3/COSMIC,  which  provided  1,500  to  2,000 
soundings  per  day,  resulting  in  GPS  RO  becoming  an 
operational  data  source  for  weather  prediction  and 
ionospheric  monitoring.  Other  missions  that  provided 
significant  quantities  of  RO  data  are  Gravity  Recovery 
and  Climate  Experiment  (GRACE- A),  METeorological 
Operations  (METOP-A),  Communications/Navigation 
Outage  Forecasting  System  (C/NOFS),  TerraS AR-X  and 
TanDEM-X.  Table  1  (taken  from  the  COSMIC  website 
[84])  shows  the  contributions  made  by  these  and  other 
missions  to  both  atmospheric  and  ionospheric  sounding. 
With  the  COSMIC  constellation  degrading,  as  it  has 
reached  its  design  life,  a  new  FORMOSAT-7/COSMIC-2 
follow-on  constellation  is  being  constructed  which  has 
multifrequency,  multi-GNSS  receivers  onboard  12  satel¬ 
lites  to  be  launched  -  six  into  low-inclination  orbits  in 
2015  and  six  into  high-inclination  orbits  in  2018  [85]. 

Far  fewer  satellites  have  been  able  to  perform  GNSS  re- 
flectometry.  The  first  was  UK-DMC  from  2003  to  2011 
[86].  Others  that  have  been  proposed  are  Techdemosat-1 
launched  July  2014  [87]  with  a  custom  reflectometry  pay- 
load;  PARIS,  a  dedicated  mission  for  GNSS-R  [35,88]; 
CYGNSS,  a  constellation  of  eight  microsatellites  [89]; 
GEROS-ISS  (GNSS  REflectometry,  Radio  Occultation  and 
Scatterometry  onboard  International  Space  Station)  [90]; 
and  Cat-2  [91].  Several  of  the  authors  are  also  from  the 
University  of  New  South  Wales  (UNSW),  which  is  prepar¬ 
ing  a  CubeSat  to  carry  both  GNSS  occultation  and  reflect¬ 
ometry  experiments  as  part  of  the  QB50  mission  [92].  As 
illustrated  in  Figure  1,  a  LEO  satellite  to  be  launched  for 
the  QB50  mission  can  receive  the  direct,  reflected  and 


Table  1  Contributions  to  atmospheric  and  ionospheric 
sounding  by  missions 


MISSION 

Total  atm  occs 

Total  ion  occs 

CHAMP 

399,968 

303,291 

CNOFS 

120,588 

0 

COSMIC 

4,039,311 

3,707,966 

GPSMET 

5,002 

0 

GPSMETSA 

4,666 

0 

GRACE 

273,013 

132,817 

METOPA 

993,084 

0 

SACC 

353,808 

0 

TSX 

276,549 

0 

Total 

6,465,989 

4,144,074 

Data  obtained  from  the  COSMIC  website  [84]  (last  updated:  Sat  Dec  7  23:25:02 
MST2013). 


refracted  GNSS  signals  to  enable  positioning,  reflectom¬ 
etry  and  occultation  applications. 

Instruments 

Even  though  primarily  designed  for  precise  navigation, 
surveying  and  geodesy,  high-end  dual-frequency  GPS  re¬ 
ceivers  have  been  used  for  more  than  two  decades  for  a 
diverse  set  of  remote  sensing  applications.  Ground- 
based  receiver  networks  such  as  that  maintained  by  the 
International  GNSS  Service  (IGS)  are  routinely  used  to 
monitor  the  tropospheric  vapour  content  and  the  iono¬ 
spheric  electron  content.  Multifrequency,  multiconstella¬ 
tion  receivers  that  are  now  available  will  further  improve 
the  spatial  and  temporal  resolution  in  the  near  future. 
With  slight  modifications  of  standard  geodetic  receivers 
(higher  data  rate,  special  observation  types),  the  so-called 
scintillation  receivers  can  be  used  to  study  the  short-term 
variability  of  the  atmosphere,  including  scintillation,  iono¬ 
spheric  storms  or  travelling  ionospheric  disturbances  [93]. 

For  satellite-borne  GNSS  applications,  various  modifi¬ 
cations  of  terrestrial  receivers  are  required,  including 
hardware  and  software  modifications  to  enhance  surviv¬ 
ability  in  a  space  environment  (temperature,  vacuum,  ra¬ 
diation)  and  to  ensure  a  proper  acquisition  of  signals 
despite  the  high  orbital  velocity  and  the  associated  Dop¬ 
pler  shifts.  On  the  other  hand,  special  hardware  and 
software  is  required  to  support  non-navigation  remote 
sensing  applications  in  space,  such  as  radio  occultation 
measurements.  Well-known  examples  include  the  Black- 
Jack  and  Integrated  GPS  Occultation  Receiver  (IGOR) 
flown  on  numerous  research  satellites  [94],  or  the  GRAS 
instrument  flown  on  MetOp  as  shown  in  Figure  6  [95]. 
These  offer  a  high-rate  sampling  (50  Hz  to  1  KHz)  of 
raw  phase  observations  or  correlator  outputs  and  open- 
loop  tracking  techniques.  A  model  of  the  expected  Dop¬ 
pler  shift  is  used  during  occultations  in  the  deep  layers 
of  the  atmosphere  to  facilitate  tracking  at  very  low  signal 
levels.  Radio  occultation  receivers  are  also  operated  with 
special  beam-forming  antennas  to  increase  the  antenna 
gain  in  directions  close  to  the  Earths  limb  [95]. 

While  the  applications  discussed  so  far  can  largely  be 
fulfilled  with  receivers  using  traditional  correlator  application- 
specific  integrated  circuits  (ASICs),  GNSS-R  calls  for  a 
higher  flexibility  in  the  receiver  design.  As  discussed  in 
the  previous  sections,  GNSS-R  may  be  used  for  a  variety 
of  remote  sensing  applications  such  as  studies  of  the  sea 
surface  height  and  roughness  or  salinity  and  soil  moisture 
from  the  properties  of  passively  reflected  GNSS  signals. 
Besides  special  antennas  (down-looking,  high-gain  or 
beam-steering  antennas  with  left-hand  polarisation)  [35], 
such  receivers  commonly  employ  a  large  number  of  slaved 
correlators  and  often  exceed  the  capabilities  of  traditional 
ASICs.  This  is  even  truer  for  new  advanced  signal  process¬ 
ing  techniques.  Recently,  two  novel  correlation  schemes 
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Figure  6  Image  of  the  GRAS  instrument  flown  on  MetOp-A.  Images  showing  the  GRAS  instrument  with  both  occultation  antennas  (front  left 

and  front  right),  the  Zenith  antenna  (front  middle)  and  the  R/F  front  ends  and  main  electronics  unit  (in  the  back). 
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have  been  introduced:  the  so-called  interferometric  pro¬ 
cessing  [35]  and  the  direct-signal-enhanced  semi-codeless 
processing  [96].  Both  methods  do  not  correlate  the 
reflected  signal  with  a  clean  replica  signal,  but  directly 
cross -correlate  the  line-of-sight  signal  received  from  a 
zenith-looking  antenna  either  with  the  reflected  signal 
from  a  nadir-looking  antenna  for  the  first  method  or  with 
a  modelled  replica  of  the  encrypted  signal  obtained  in  a 
semi-codeless  way  for  the  second  method.  Key  advantages 
of  these  new  schemes  are  that  no  prior  knowledge  of  the 
incoming  signal  is  required  and  there  is  a  significant  im¬ 
provement  in  the  measurement  precision  and  the  carrier 
to  noise  ratio.  For  altimetric  applications  in  particular, 
which  are  probably  the  most  demanding  applications  in 
terms  of  the  required  instrument  hardware  and  processing 
capabilities,  this  is  of  crucial  benefit. 

The  disadvantage  of  these  novel  methods  is  that  the  cor¬ 
relation  requires  a  substantially  increased  signal  process¬ 
ing  effort.  To  cope  with  these  needs,  a  software-defined 
radio  approach  is  widely  applied.  In  the  extreme  case,  pure 
software  receivers  have  been  employed  by  various  research 
groups  for  offline  processing  of  the  pre-recorded  radio  fre¬ 
quency  (RF)  signals,  while  field-programmable  gate  arrays 
(FPGA)  or  fast  digital  signal  processors  (DSPs)  are  com¬ 
monly  used  in  real-time-capable  reflectometry  receivers. 
Examples  of  the  latter  type  of  receivers  include  the  GOLD- 
RTR  [97]  and  SSTLs  SGR-ReSI  [98]  instruments  developed 
by  IEEC,  Barcelona,  and  SSTL,  Surrey,  respectively. 

Challenging  issues  and  future  directions 

GNSS  has  a  global  space-based  positioning,  navigation 
and  timing  (PNT)  capability.  There  are  many  techno¬ 
logical  challenges  for  GNSS  remote  sensing  to  come  close 
to  matching  the  extraordinary  success  of  GNSS-PNT, 


especially  in  the  case  of  the  GNSS-R  technique.  Satellite- 
borne  GNSS  remote  sensing  must  demonstrate  its  value 
as  a  reliable,  continuous,  high-quality  sensing  technology. 
With  the  launch  of  multiple  LEO  satellites  with  GNSS- 
RO-  and  GNSS-R-capable  instrumentation,  it  seems  we 
are  at  last  close  to  a  renaissance.  But  are  more  LEO  mis¬ 
sions  and  better  algorithms  for  geophysical  parameter  ex¬ 
traction  sufficient?  In  the  case  of  GNSS-RO,  the  answer  is 
yes  because  the  meteorological/climate  community  now 
assimilates  GNSS-RO  products  into  their  operational  and 
research  systems  and  is  eager  for  more  [11-15].  However, 
GNSS-R  is  still  very  much  a  novelty  technology  as  far  as 
the  geosciences  community  is  concerned  because  current 
results  are  generated  by  one  or  more  isolated  academic 
group  and  not  organised  into  independently  verifiable 
products.  Nevertheless,  as  mentioned  earlier,  it  is  a  fact 
that  the  Plate  Boundary  network  is  being  used  in  the  USA 
for  a  range  of  GNSS-R  applications. 

The  International  Association  of  Geodesy  (IAG)  has 
been  very  successful  in  launching  technique-specific  ser¬ 
vices  (see  [99]).  These  services  aid  other  geoscientists  in 
their  research,  as  well  as  supporting  important  applica¬ 
tions  in  the  wider  community.  Examples  of  the  latter  are 
the  contribution  of  the  International  GNSS  Service  (IGS) 
to  precise  positioning  [100]  and  the  International  Earth 
Rotation  and  Reference  Systems  Service  (IERS)  to  the 
International  Terrestrial  Reference  Frame  (ITRF)  [101]. 
The  IGS  does  generate  troposphere  and  ionosphere  par¬ 
ameter  products  [102];  however,  they  are  based  on  ob¬ 
servations  made  from  the  global  ground-based  GNSS 
tracking  network  (Figure  7),  not  from  satellite  plat¬ 
forms.  There  are  no  geodetic  services  producing  GNSS 
remote  sensing  products  on  a  continuous,  synoptic 
basis.  From  the  IAG  s  perspective,  the  challenging  issues 
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Figure  7  IGS  observations  made  from  the  global  ground-based  GNSS  tracking  network.  One  of  the  reasons  for  the  success  of  IAG  services 
are  the  permanent  geodetic  observing  networks  such  as  IGS's  global  GNSS  tracking  network,  providing  accurate,  reliable  data  products  for 
science  and  society. 
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are  related  to  operationalising’  GNSS  remote  sensing,  in 
all  its  forms,  so  that  the  remotely  sensed  geophysical  pa¬ 
rameters  (atmospheric,  oceanic,  wind,  soil  moisture,  bio¬ 
mass,  etc.)  are  included  in  the  suite  of  geodetic  outputs  of 
the  Global  Geodetic  Observing  System  (GGOS)  [103]. 

The  non-technical  challenges  therefore  include  the  evo¬ 
lution  of  current  GNSS  remote  sensing  science  missions 
to  operational  services  by  increasing  the  number  of  LEO 
satellites  equipped  with  GNSS  receivers  to  satisfy  coverage 
requirements  in  space  and  time;  standardising  data  for¬ 
mats,  instrumentation,  and  calibration,  pre-processing  and 
analysis  procedures;  and  establishing  a  coordinating  agency 
or  authority  so  as  to  ensure  continuous,  high-quality  prod¬ 
uct  generation  and  dissemination.  The  IAG  has  consider¬ 
able  experience  with  geometry  technique-based  services 
(such  as  the  IGS,  ILRS,  IVS  and  IDS)  [104];  however,  it  has 
not  yet  established  any  ‘geodetic  imaging’  services  based 
on  technologies  such  as  synthetic  aperture  radar  (SAR), 
satellite  radar  altimetry,  lidar  or  GNSS  remote  sensing 
(GNSS-RO,  GNSS-R).  One  of  the  future  challenges  is  to 
address  this  shortcoming,  so  that  GNSS  remote  sensing 
can  be  recognised  as  a  geodetic  technique  that  is  making 
critical  contributions  to  science  and  society. 

Conclusions 

This  paper  provided  an  overview  of  GNSS  remote  sensing 
(GNSS-RO  and  GNSS-R)  with  a  focus  on  six  different  is¬ 
sues.  Atmospheric  sensing  using  GNSS-RO  measurements 
has  been  extensively  investigated  for  determining  atmos¬ 
pheric  parameters  such  as  water  vapour,  temperature,  pres¬ 
sure  and  humidity.  Significant  benefits  have  been  shown 
for  the  southern  hemisphere  where  there  is  limited  access 
to  ground-based  observations,  and  GNSS-RO  has  been 


operationally  used  by  major  international  weather  fore¬ 
casting  models.  GNSS-R,  on  the  other  hand,  which  has 
recently  gained  more  attention  from  research  organisa¬ 
tions  and  universities,  can  be  used  to  perform  remote 
sensing  tasks  related  to  ocean  observation  and  land  ap¬ 
plications,  retrieving  a  wide  range  of  geophysical  param¬ 
eters  including  sea  surface  wind  speed  and  direction, 
sea  surface  height  and  roughness,  soil  moisture,  bio¬ 
mass,  snow/ice  depth,  coverage  and  density.  The  most 
important  past  and  planned  satellite  missions  primarily 
associated  with  GNSS-RO  and  GNSS-R  applications 
were  summarised.  More  GNSS-RO  missions  have  been 
launched,  while  a  few  GNSS-R  missions  are  planned. 
The  state-of-the-art  characteristics  of  the  GNSS  receiver 
were  discussed,  which  is  the  key  equipment  for  GNSS 
remote  sensing.  Finally,  a  number  of  challenging  issues 
related  to  GNSS  remote  sensing  products  and  services 
were  highlighted. 
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